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ABSTRACT 

We report on very-long-baseline-interferometry (VLBI) monitoring observations of the low-luminosity active galactic nucleus 
(LLAGN) in the galaxy M81 at the frequencies of 1.7, 2.3, 5.0, and 8.4GHz. The observations reported here are phase-referenced to 
the supernova SN 1993J (located in the same galaxy) and cover from late 1993 to late 2005. The large amount of available observa- 
tions allows us to study the stability of the AGN position in the frame of its host galaxy at different frequencies and chromatic effects 
in the jet morphology, together with their time evolution. The source consists at all frequencies of a slightly resolved core and a small 
jet extension towards the northeast direction (position angle of ~65 degrees) in agreement with previous publications. We find that the 
position of the intensity peak in the images at 8.4 GHz is very stable in the galactic frame of M 8 1 (proper motion upper limit about 
lOyuas per year). We confirm previous reports that the peaks at all frequencies are systematically shifted among them, possibly due to 
opacity effects in the jet as predicted by the standard relativistic jet model. We use this model, under plausible assumptions, to estimate 
the magnetic field in the jet close to the jet base and the mass of the central black hole. We obtain a black-hole mass of ~2 x 10 7 M Q , 
comparable to estimates previously reported using different approaches, but the magnetic fields obtained are 10 3 - 10 4 times lower 
than previous estimates. We find that the positions of the cores at 1.7, 2.3, and 5.0 GHz are less stable than that at 8.4 GHz and evolve 
systematically, shifting southward at a rate of several tens of yuas per year. The evolution in the jet orientation seems to be related to 
changes in the inclination of the cores at all frequencies. These results can be interpreted as due to a precessing jet. The evolving jet 
orientation also seems to be related to a flare in the peak flux densities at 5.0 and 8.4 GHz, which lasts ~4 years (from mid 1997 to 
mid 2001). An increase in the accretion rate of the black hole, and its correlation with the jet luminosity via the disk-jet connection 
model, seems insufficient to explain this long flare and the simultaneous evolution in the jet orientation. A continued monitoring of 
the flux density and the jet structure evolution in this LLAGN will be necessary to further confirm our jet precession model. 
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The galaxy M81 (NGC3031), located at a distance of 3.63 + 
0.34 Mpc (Freedman et al. I1994K hosts the next closest ac- 
tive galactic nucleus (AGN) after Centaurus A. Its high dec- 
lination (J2000.0 coordinates a = 09 h 55 m 33.173 s and 5 = 
69° 03' 55.062") also makes M81 a perfect target for global 
_._ H VLBI observations using the most sensitive radio telescopes in 
• • . the Northern Hemisphere. 

The flux density of the AGN in M 8 1 (hereafter, M 8 1 *), has 
h>^ relatively large variability, from X-rays (Ishisaki et al. 119961 
$_i ' Page et al. 1994) to radio wavelengths (Ho et al. 119991 and 
Bietenholz, Bartel & Rupen 120001 at centimeter wavelengths; 
Schodel et al. 120071 at millimeter wavelengths). The variabil- 
ity time scale runs from years to days, indicating a contribution 
from a very compact source to the emission at all frequencies. 
Indeed, even from VLBI observations at 43 GHz (Ros & Perez- 
Torres 2008) and space-VLBI observations at 5 GHz (Bartel & 
Bietenholz 120001) the emitting structure could not be fully re- 
solved. The broadband luminosity of M 8 1 * is relatively weak 
(see, e.g., Ho et al. 1999) from radio (around 100 mJy, i.e., on the 
order of 10 37 ergs~') to the optical and X-rays (~ 10 40 ergs _1 ). 
Therefore, it is classified as a low-luminosity AGN (LLAGN). 

The accretion mechanisms in LLAGNs are not clear. Since 
relativistic X-ray lines are not clearly detected in these objects 
(e.g., Reynolds et al. 120091 1. radial flows or advection mecha- 



nisms might dominate the accretion, in contrast to the more lu- 
minuos AGN. In any case, if the X-ray emission traces the in- 
nermost inflow from the disk to the central engine and the radio 
emission traces the jet power, there should be a relationship be- 
tween luminosities in both X-rays and radio for all AGNs (i.e., 
including LLAGNs), which should also be related to the accre- 
tion rate. This is the basis of the fundamental-plane model of 
black-hole accretion (e.g., Merloni, Heinz, & di Matteo 2003). 
However, the parametric space of very-low accretion flows in 
AGNs (e.g., LLAGN) is not fully characterized observationally. 
Although Markoff et al. (2008 ) and Miller et al. (2010) were able 
to fit average broadband spectra of M 8 1 * in the frame of the 
fundamental-plane model of disk-jet connection, no clear corre- 
lation between the variabilities in the X-rays and radio was ob- 
served unlike in microquasars (Mirabel et al. 1998) and AGNs 
(Marscher et al. 120021 ). In the present paper, we report on a very 
long flare of M 8 1 * at radio frequencies, beginning around mid 
1997, which lasted nearly four years. As shown in this paper, an 
intriguing result related to this flare is that it seems to be related 
to the evolution in the jet orientation, thus making it difficult ex- 
plain in the frame of the disk-jet connection between accretion 
flow and jet power. 

At radio frequencies, M81* has a slightly inverted spec- 
trum with an spectral index a - +0.3 (S oc v +a ), up to a 
turnover frequency of ~200GHz (Reuter & Lesch [l996l ). If due 
to synchrotron self-absorption in the jet, such a high turnover 
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frequency, suggests a very dense plasma and a high magnetic 
field, and/or a rather flat energy distribution of the relativis- 
tic synchrotron-emitting particles (see, e.g., Sect. 4.1 of Reuter 
& Lesch 1996). In addition, it is worth noticing that the spec- 
tral behavior of M 8 1 * is similar to that of the central engine 
in our Galaxy, Sgr A* (Duschl & Lesch |1994| l, thus suggesting 
that similar physical mechanisms should be present in the ac- 
celeration of the relativistic particles and their emission at radio 
frequencies. In addition, the mass of its associated black hole 
(Tl x 10 7 M Q , De vereux et al. [20031 5.5+ 3 6 x 10 7 M o , Schorr 
Miiller et al. l201U makes M 81* an interesting intermediate ob- 
ject between Sgr A* and the more powerful AGN. 

The explosion of the radio-luminous supernova SN 1993 J 
in M81, which happend around 28 March 1993 (e.g., Weiler 
et al. 120071 1. triggered an intense campaign of VLBI observa- 
tions in which M 8 1 * was selected as the phase calibrator of the 
SN 1993J visibilities. These observations allow us to study the 
evolution in the structure and position of M 8 1 *, in the frame of 
its host galaxy and at different frequencies and times. Results for 
some of the M81* observations at 8.4 GHz between years 1993 
and 1997 are reported in Bietenholz, Battel & Rupen (2000), and 
results for the astrometry analysis of a subset of observations at 
all frequencies between years 1993 and 2002 were reported in 
Bietenholz et al. ([2004). These authors conclude that the proper 
motion of M81*, relative to the shell center of SN 1993J, was 
consistent with zero up to the precision level of the observa- 
tions ( 1 1 .4 ± 9.3/iasyr~ 1 ). The contributions of galactic rotation 
and/or any peculiar motion of the SN 1993 J progenitor were also 
discarded to contribute to a detectable proper motion of M 8 1 * 
relative to the supernova. These authors also report a frequency- 
dependent shift of the peak of brightness of M81*, which they 
interpret as opacity effects in the jet of the AGN, in agreement 
with the standard relativistic jet model (Blandford & Konigl 
[19791 . 

In this paper, we report on the analysis of all the VLBI obser- 
vations of M81* taken in the phase-referencing observing cam- 
paign of SN 1993 J, from late 1993 to late 2005. In Sect. [2] we de- 
scribe our observations. In Sects.[3]and[4]we report the main re- 
sults obtained (source morphology and multi-frequency astrom- 
etry). In Sect. [5] we analyze the astrometry results in terms of 
the standard jet interaction model (Blandford & Konigl [T979I ), 
and in the last section we discuss the detection of a changing jet 
orientation, which we interpret as due to jet precession. 



2. Observations and data reduction 

We analyzed the VLBI data from the observing campaigns 
led by N. Bartel (York University, Canada) and J.M. Marcaide 
(University of Valencia, Spain), from late 1993 to late 2005 (see 
Table 1 of Marti- Vidal et al. 120 11 al and Table 1 of this paper). 
SN 1993 J was observed using M81* as the phase calibrator at 
the frequencies of 1.7, 2.3, 5.0, 8.4, 15, and 22 GHz. (There are 
only a few low-quality observations during years 1993 and 1994 
at the two highest frequencies, which have not been used in our 
analysis.) The results of this observing campaign of SN 1993J 
are reported in several publications (e.g., Bartel et al. |2002, 
Marcaide et al. [20091 Marti- Vidal et al. 1201 la!201 lbt and ref- 
erences therein), where the technical details of all these obser- 
vations are given. Basically, one epoch typically lasted ~12h, 
with the participation of 10-15 antennas. The whole Very Long 
Baseline Array (VLBA), the Very Large Array (VLA) in phased- 
array mode, and some stations of the European VLBI Network 
(EVN) took part in most of the observations. The recording rate 



was set to either 128 or 256 Mbit s , depending on the technical 
capabilities of the recording systems at each epoch. 

To obtain the calibrated visibilities of M 8 1 * at each epoch, 
we proceeded as described in Marti- Vidal et al. J201 lal >. The 
Astronomical Image Processing System (aips) of the National 
Radio Astronomy Observatory (NRAO) was used to calibrate 
the data. Fringe-finder scans were used to align the phases be- 
tween the different observing sub-bands and, afterwards, global 
fringe fitting (Alef & Porcas 1986) was applied to the M 81* vis- 
ibilities. The resulting gains were interpolated to the SN 1993J 
scans. The amplitude calibration was based on gain curves and 
system temperatures measured at each station, with the excep- 
tion of the phased VLA (where estimates of the source flux den- 
sities are necessary, since the amplitude calibration is based on 
measurements of the ratio of effective antenna temperature to 
system temperature). Then, a refinement of the amplitude gain 
estimates was performed via self-calibration of the M 8 1 * visi- 
bilities, using the aips task CALIB. 



3. Results for source morphology 

We show in Fig. Q] some representative examples of the VLBI 
images obtained at all the observing frequencies. The images of 
M 8 1 * consist of a compact source, which dominates the emis- 
sion, and a weak extension towards the northeast. These images 
agree with the results reported in Bietenholz, Bartel, & Rupen 
(2000). 



3.1. Model fitting 

For the analysis of evolution in the M 8 1 * structure, we char- 
acterized the model emission at each epoch in a similar way to 
what was modeled in Bietenholz, Bartel, & Rupen (2000). We 
used an elliptical Gaussian to model the central dominant emis- 
sion and a point source to model the weak jet extension. For 
those epochs in which the point-source model could not describe 
the M 81* jet extension satisfactorily, we used a circularly sym- 
metric Gaussian. The point source used to succesfully model the 
jet extension at some epochs is indicative of a compact emission 
away from the core (i.e., a hot spot in the jet), at least up to the 
sensitivity limit achieved in the images. 

The defining parameters of the geometrical model described 
in the previous paragraph were estimated for each epoch by 
model fitting to the visibilities, as implemented in the program 
difmap. We show the best-fit model parameters at all frequencies 
and epochs in Table 1 and Figs. EJ [3] EJ and [5] For the model 
of the core emission, the parameters are peak intensity, position 
angle (both shown in Fig. 0, major axis, and axis ratio (both 
shown in Fig. 0. For the model of the jet extension, the pa- 
rameteres are separation from the peak of the core, orientation 
with respect to the core (both shown in Fig. @}, and flux den- 
sity (shown in Fig. |5J, We notice that the best-fit parameters of 
the core Gaussians are insensitive to the presence of the addi- 
tional model component used to describe the jet extension, due 
to the much smaller contribution of that component to the over- 
all structure of the source. (We notice that the flux density of the 
jet extension is typically below 7% of the peak flux density of 
the core, see Fig. [5]) We have discarded some outliers in Table 
1 to generate Figs. [2] [3] 13 and[5] Basically, we discarded those 
fits in which the final axis ratio of the core Gaussian was zero 
(i.e., the peak emission could not be satisfactorily described as 
a 2-dimensional Gaussian) and the few epochs in year 1993 that 
result in negative position angles. 
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Fig. 1. Images of M 8 1 * obtained from the VLBI epochs at (a) 1 .6 GHz, (b) 2.3 GHz, (c) 5.0 GHz, and (d) 8.4 GHz observed around 
mid November 2000 (see Table 1). Contours are located at 3, 3 V3, 9, ... of the rms of the images. The FWHM of the convolving 
beam is shown in the bottom left corner of each image. 



The typical values for the uncertainties in the model param- 
eters of the core Gaussians were estimated using the following 
approach. For some representative epochs, we found the peak in 
the image of the residuals corresponding to the best-fit model. 
(We call this quantity <x.) Then, we generated a synthetic im- 
age of the best-fit model, free of noise from the visibilities. 
We call this image /best- Then, to check the range of parame- 
ter values that still lead to satisfactory fits, we changed the val- 
ues of the parameters (one parameter at a time), generated the 
corresponding noise-free model images (we call them /,•), and 
computed the maximum deviation between 4 es t an d a ^ h (i- e -' 



cri - max (|/best - Ii\))- The parameter uncertainties were esti- 
mated as the changes in the parameter values such that a-j — 3cr. 
Following this approach, the resulting uncertainties are those 
that translate into changes in the model that are three times 
greater than the highest noise contribution to the image. The 
typical uncertainties estimated in this way are given in Table [2] 
These are the uncertainties assigned to the parameters at each 
epoch. However, the uncertainty in the average of the parameter 
values for a given set of epochs is set to be equal to the standard 
deviation of the parameter values, to reflect the natural scatter in 
the parameters from epoch to epoch. 
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Fig. 2. Peak flux density (left) and position angle (right) of the elliptical Gaussians used to model the core emission of M 8 1 * at all 
frequencies (see Table 1 and text). Typical uncertainties are given in Table|2] 
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Fig. 3. Major axis (left) and axis ratio (right), i.e., minor axis in units of major axis, of the elliptical Gaussians used to model the 
core emission of M 81* at all frequencies (see Table 1 and text). Typical uncertainties are given in Table|2] 
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Fig. 4. Separation from the core (left) and orientation relative to the core (right) of the jet extension of M 8 1 * at all frequencies. 
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Table 2. Typical uncertainties in the model parameters of the 
core Gaussians. 
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Fig. 5. Flux density of the M81* jet extension relative to the 
peak intensity of the core. 



terpret the jet extension as an optically-thin region in the jet at 
all frequencies. In such a case, the distance between the central 
AGN engine and the (optically-thin) jet extension should be the 
same at all frequencies. The distance between the core and the jet 
extension should therefore systematically decrease at lower fre- 
quencies if the distance between the central AGN engine and the 
core increases (as we show in Sect. 15.21 ). The strange behavior 
seen in Fig. [4] (left) may indicate that the jet extensions at dif- 
ferent frequencies are spatially located in different regions along 
the jet, although other problems related to the very different spa- 
tial resolutions achieved at different frequencies (i.e., blending 
of jet components at lower resolutions) could also be affecting 
the results of our model fitting. 

Indeed, we notice that the reported jet extensions at 8.4 and 
5 GHz are located within the core at 1 .7 GHz (the typical core 
size at 1.7 GHz is ~2-2.5 mas). Therefore, the jet extension at 1.7 
GHz (which is located at a distance of ~3-4 mas from the core 
at the same frequency) must be a physically distinct component, 
indicating that the successive jet extensions at all frequencies are 
indeed tracing different parts of a continuous jet. Figure|4](right) 
shows that the position angle of the jet extension is similar at 
all frequencies, thus suggesting a rather straight jet in the region 
of the emission at all our observing frequencies. This result is 
supported in the following section, but somewhat in conflict with 
the results reported in Sect. 15.31 

Finally, Fig. [5] shows the flux density of the jet extension, 
in units of the core flux density, for all epochs and frequencies. 
The typical flux density of the extension is lower than ~7% of 
the core flux density, with the exception of only two epochs: 
an epoch at the end of year 1996 (around 10% of the core flux 
density at 8.4 GHz) and another at the end of year 2005 (around 
18% of the core flux density at 5.0 GHz). 



Figure [2] (left) shows how the peak intensity at 5.0 and 
8.4 GHz of the fitted Gaussians increases from June 1997 to July 
1998 and then begins to decrease through June 2001. The max- 
imum peak intensity is a factor ~2 of the typical peak intensi- 
ties outside this four-year-long flare. Radio-variabilities of sev- 
eral years have been reported for several LLAGN and Seyfert 
galax ies (Falcke et al. [200T1 Nagar et al. 120021 Mundell et al. 
2009). However, in this case, as shown below, we also notice 
hints of a simultaneous evolution in the orientation of the fitted 
core-Gaussians, which increase the position angle of their major 
axes from ~60 to ~70 degrees (Fig. [2] right) during the flare. In 
Sect. 15.31 we also discuss the possible correlation between this 
flare (together with the changes in the orientation of the fitted 
core-Gaussians) and our astrometry results at the different fre- 
quencies. 

Figure [3] shows that the size of the fitted core-Gaussians 
(i.e., basically the size of the core in the images) is system- 
atically larger at lower frequencies (varies from ~0.5mas at 
8.4 GHz to ~2mas at 1.7 GHz). The ratio between minor and 
major Gaussian axes is also higher at lower frequencies, thus in- 
dicating that the core regions at lower frequencies are relatively 
broader. A larger core size at lower frequencies, combined with 
the frequency-dependent core shift reported in the following sec- 
tions, indicates an opening angle in the jet, as discussed in Sect. 
El 

Figure [4] (left) shows how the separation between the peak 
intensity and the jet extension increases at lower frequencies, 
with the exception of the data at 2.3 GHz, which gives a larger 
separation than at 1 .7 GHz. This result is intriguing if we in- 



3.2. Extended emission at 1 .7 GHz 

In the most sensitive low-frequency VLBI observations (those at 
1.7 GHz with large arrays and long observing times), we detect 
the signature of an extended emission centered at ~ 15 mas from 
the core with a peak intensity of ~0.7 mjy beam -1 . The position 
angle of the peak of this extended structure is ~ 65 ° (north- 
east direction). We show in Fig. [6] an image with this extended 
emission. We notice that the intensity of this component is so 
low that natural weighting of the visibilities^, plus an additional 
downweight of the longest baselines (a Gaussian taper with half- 
width at half-maximum, HWHM, of 30 MA or lower) is neces- 
sary to detect the extended emission in the images. Since it is not 
detected at the higher frequencies, we can only set a lower ab- 
solute value for the spectral index of this extended component, 
based on the rms of the most sensitive observations at the higher 
frequencies. We find for this extended component that a < -0.78 
(being the flux density S oc v +a ). 

The similarity of position angles between the jet components 
reported in the previous section and this more extended emission 
points toward a rather straight jet in M 8 1 * or a jet with no strong 
projected bendings up to a distance of at least ~ 15 mas (about 
0.26 pc of projected linear scale). 



1 Each pixel in the Fourier plane is weighted according to the scatter 
of the visibilities inside it. This weighting scheme improves the sensi- 
tivity, but decreases the resolution. 
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4. Results for multifrequency astrometry 

4. 1 . Frequency-dependent core shifts 

In Fig. [7j we show the position of the peak of brightness of 
M81* relative to the J2000.0 coordinates a = 09 h 55 m 33.173 s 
and 6 = 69° 3' 55.062", and with the center of the SN 1993J ra- 
dio shell as a position reference. These coordinates correspond 
to the nominal position of M 8 1 * used in the correlation of the 
data at 5 GHz taken on 30 May 1998. We set the position of the 
M 8 1 * core at this epoch as our coordinate origin. Therefore, all 
the other phase-referenced observations correlated using differ- 
ent nominal positions for SN 1993J and/or M81* were accord- 
ingly shifted. In this figure, we show the results for the obser- 
vations at 8.4, 5.0, 2.3, and 1.7 GHz, which are also given in 
Table 1 . The uncertainties are estimated as the half width at half 
maximum (HWHM) of the interferometric beam, divided by the 
dynamic range of the image of SN 1993 J corresponding to each 
epoch. The uncertainties computed in this way are related to the 
statistical contribution of the image noise in the location of the 
SN 1993J radio shell. They are also a good estimate of the as- 
trometric precision of compact sources observed with the phase- 
referencing technique (e.g., Marti- Vidal et al. 1201 Oi l. 

Figure [7j shows the dependence of the position of the peak 
of brightness of M 8 1 * on the observing frequency, taking the 
center of the SN 1993 J shell (at each frequency) as a position 
reference. We can use this figure to study the shift of the peak 
of brightness of M 8 1 * as a function of frequency in the con- 
text of the relativistic jet model (e.g., Blandford & Konigl 1979, 
Lobanov 1998). From the time average of the M81* positions 
shown in Fig. [7] at each frequency, we derived the relative shifts 
shown in Table [3] (Row 1) for all pairs of frequencies. 
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Fig. 7. Shifts in right ascension and declination of the peak in- 
tensity of M 8 1 *, relative to a fixed position on the sky (see text) 
and using the center of the SN 1993 J radio shell as a position 
reference. 



There may be, however, some systematic effects in the data 
that should be taken into account, or be observationally dis- 
carded, before we take the shifts in Table[3](Row 1) as true shifts 
of the brightness peak of M 8 1 *. On the one hand, the positions 
shown in Fig. [7] are obtained from fits of a circularly symmet- 
ric shell to the SN 1993J visibilities. The shell of SN 1993J does 
not have a perfect circular symmetry; there are deviations with 
respect to a perfect circular source in both the shape of the shell 
(of the order of 2%) and the intensity distribution inside it (of the 
order of 20%) (Bart el et al . 120021 Bietenhol z, Bart el & Rupen 
2000; Marcaide et al. 120091 Marti- Vidal et al. 120 Hat . Moreover, 
the inhomogeneities in the shell evolve systematically in time. 
Therefore, and especially for the latest epochs when the super- 
nova is so extended and weak, there could be unexpected sys- 
tematics in the estimate of the shell center. On the other hand, 
although the proper motion of SN 1993 J in the galactic coordi- 
nates might be small (due to the peculiar motion of the supernova 
progenitor), its integrated effect during more than a decade could 
still leave a small fingerprint in the position estimates of M 8 1 * 
(see the next subsection). 

We decided to apply a different approach to study the 
frequency-dependent shift of the M 8 1 * brightness peak, which 
solves, in principle, the drawbacks of using the supernova shell 
center at each epoch as a position reference. For the subset 
of epochs where there are quasi-simultaneous observations of 
the supernova at more than one frequency, we convolved the 
CLEAN models of SN 1993 J at each epoch with beams of the 
same size, namely, 0.5 times the shell radius at 5 GHz reported 
in Marti- Vidal et al. (120 1 1 at . We then cross-correlated the re- 
sulting CLEAN supernova images at different frequencies and, 
from the peak of the cross-correlation image, we estimated the 
shift of M 81* between each pair of frequencies. In Fig. [8j we 
show two images of SN 1993J (at 5 and 8.4 GHz), correspond- 
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Table 3. Shifts in the M 81* brightness peak for several pairs of observing frequencies. 



Method 

Shell center" 
Cross-corr.* 



1.7-5.0 



1.7-8.4 



Frequency pair (GHz) 

5.0-8.4 1.7-2.3 



2.3-5.0 



2.3-8.4 



0.57: 
0.68: 



0.21 
0.22 



0.74: 
0.77: 



0.19 
0.50 



0.18: 
0.21: 



Shift (mas) 

0.13 0.31: 
0.10 0.35: 



0.17 
0.05 



0.42: 
0.43: 



0.14 
0.25 



0.59: 
0.65: 



0.09 
0.14 



(fl) Taking the estimated center of the SN 1993J radio shell (at each epoch and frequency) as position reference. {b) From the shift in the 
cross-correlation of SN 1993J images at different frequencies and similar epochs. 




Relative Dec (mas) 



Relative Dec (mas) 



Relative Dec (mas) 



Fig. 8. Left and center: CLEAN images of SN 1993J observed at 5.0GHz and 8.4 GHz, respectively, around September 1996. The 
color scale runs from black (minimum intensity) to white (maximum intensity). Both images are phase-referenced to M 8 1 *. Right: 
image cross-correlation (black=-0.02; whiter 1.00). Notice the small shift in the correlation peak from the image origin. 



ing to epochs around September 1996, together with their image 
cross-correlation. The small shift in the correlation peak from 
the center of the image can be appreciated (~ 0.2 mas roughly 
to the west). In Table [3] (Row 2), we present the average shifts 
of the M 8 1 * brightness peak estimated using this approach. The 
uncertainties are computed as the standard deviation of the aver- 
ages, to reflect the natural scatter of the data. For all cases, the 
error bars are larger than those determined from model fitting. 
We used a maximum gap of 15 days between two epochs to con- 
sider them as quasi-simultaneoufl and also required a minimum 
dynamic range (i.e., peak intensity over root-mean-square of the 
residual images) of ten in both SN 1993J and M 81*, to use the 
data of the corresponding epoch. We notice that the shifts com- 
puted from the shell-center estimate (Row 1 in Table [3]), which 
are based on model fitting to the visibilities, are very similar to 
those computed from the cross-correlation approach (Row 2 of 
the table), which is a method completely different from model 
fitting. Both independent methods therefore give consistent re- 
sults. 

We show in Fig. [9] the resulting shifts of the M81* bright- 
ness peak for several pairs of observing frequencies, obtained 
with the cross-correlation approach. Any relative motion be- 
tween SN 1993 J and M81* during the whole observing cam- 
paign, and any systematics coming from the fit of a circular 
shell to the evolving supernova structure should not affect the 
results shown in this figure. From this figure, there is no clear, 
unambiguous, systematic time evolution in the shifts during the 
time period covered by our observations, save for the pair of fre- 



quencies 1 .7-5 GHz at late epochs, with a decrease in relative 
declination and a slight increase in relative right ascension. This 
evolution could be indicative of precession of the intensity peak 
at 5 GHz with respect to the peak at 1 .7 GHz. However, these lat- 
est astrometry results depend on the images of SN 1993 J at late 
epochs, which are the noisiest. The evidence of jet precession re- 
sulting only from Fig. [9] should, thus, be approached with care. 
However, in Sect. I5.3l we provide evidence of jet precession re- 
sulting from the combination of the multifrequency astrometry 
and the morphological evolution of the jet. 

4.2. Position stability ofM81* 

Figure [7] shows a frequency-dependent shift in the brightness 
peak of M 8 1 *. In Table|U we report the proper motions of those 
peaks in right ascension, /u a , and declination, fig, considering all 
the available data at each frequency. The position of the inten- 
sity peak at 8.4 GHz is the most stable one, with a proper motion 
consistent with zero at a level of less than 10 yt/as year 1 . In the 
other cases but 1.7 GHz, the proper motion fitted in declination 
is negative and above 3<x. 

The proper motion in right ascension is consistent with zero 
at all frequencies but 5.0GHz, where a decreasing right ascen- 
sion of lOyuasyear -1 is found at a 3cr level. These results dif- 
fer from those reported in Marti- Vidal et al. (1201 la) , based on 
the same astrometric dataj, but converted into separations and 
position angles (measured north through east). In that publica- 
tion, we reported that no change in the separation between the 



2 From the model of Marti- Vidal et al. (201 la), the expansion of the 
supernova during 15 days is only ~ 30/jas on day 500, and less at later 
times, which is well below the VLBI resolution for all observation fre- 
quencies reported here. 



3 We notice, however, that an uncertainty of 0. 1 mas was assigned to 
all the position estimates of the peak intensity of M 8 1 * in that publica- 
tion, which differs from the uncertainties used in the present paper. 
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Fig. 9. Time evolution of the shift in the intensity peak of M 8 1 * 
observed at different pairs of frequencies (peak at the lowest fre- 
quency referred to the peak at the highest frequency). 

Table 4. Proper motion of the intensity peak of M 8 1 *. 



Freq. (GHz) 


fi a (mas year ' ) 


/i,5 (mas year ' ) 


1.7 


0.016 ±0.018 


-0.010 ±0.018 


2.3 


0.008 ± 0.022 


-0.076 ± 0.022 


5.0 


-0.010 ±0.003 


-0.034 ± 0.003 


8.4 


0.000 ± 0.007 


-0.007 ± 0.007 



intensity peak and our coordinate origin was detected at any fre- 
quency, with the exception of the data at 5 GHz beyond July 
2001 (i.e., -3000 days after the explosion of SN 1993J). This 
is, indeed, a correct result, but a change in separation consis- 
tent with zero was then interpreted as a proper motion that is 
also consistent with zero. Such an interpretation does not con- 
sider the contribution of any possible rotation with respect to 
the reference point, which may contribute, however, to the right- 
ascension and declination proper motions shown in Table [4] In 
either case, the small proper motions reported in the present pa- 
per translate, at most, into shifts with respect to the SN 1993 J 
shell center of only 1 - 2% of the shell size at any epoch. These 
shifts are, indeed, much smaller than the scatter in the positions 
of the M 8 1 * peak, which is similar to the shell-size uncertainties 
reported in Marti- Vidal et al. (2011a) (see a detailed discussion 
on this matter, and its potential impact in the reported expansion 
curve of SN 1993J, in the last paragraph of their Sect. 3.1). 

As a complementary plot, we show in Fig. [10] the evolution 
on the sky plane of the position of the brightness peak at all fre- 
quencies and epochs. On the one hand, different frequencies are 
plotted using different symbols; on the other hand, and for rea- 
sons of clarity, the results at different epochs are plotted using 
different colors. Although there are some outsiders and the scat- 
ter in the data is large, Fig.[10]shows a clear time evolution in the 
overall orientation of the jet, having smaller position angles (i.e., 
orientation towards north) at earlier epochs (i.e., reddish colors) 
compared to those at later epochs (i.e., greenish and blueish col- 
ors). The data at 8.4 GHz are only present at early epochs (until 
November 2000), while data at 1.7 GHz are only present at late 
epochs (beginning on November 1998). 
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Fig. 10. Position of the brightness peak of M 8 1 * at different fre- 
quencies (i.e., different symbols) and different epochs (i.e., dif- 
ferent colors). The points plotted are taken from columns Aa and 
A<5 of Table 1. 



5. Discussion 

5. 1 . Spectral -index images 

The radio emission from SN 1993 J at all frequencies comes 
from the very same location: the shocked circumstellar region. 
Therefore, the use of SN 1993 J as a position reference, either 
using the estimated shell center or cross-correlating the images 
at different frequencies, allows us to build spectral-index images 
of M81* free of assuming any co-spatial optically-thin jet com- 
ponent at different frequencies, which is, by far, the most com- 
monly used assumption in the study of chromatic effects in AGN 
jets (e.g., Kovalev et al. 120081 1. 

We used the shifts in the SN 1993 J image cross-correlations 
(see Sect. 14. It to align the M 8 1 * images at different epochs and 
for different pairs of frequencies. The results for two representa- 
tive epochs are shown in Fig. QT| In each image, we convolved 
the CLEAN models at both frequencies using the beam corre- 
sponding to the lowest frequency. To avoid spurious contribu- 
tions to the spectral-index maps, we clipped those points in the 
image with an intensity lower than 2% of the peak intensity at 
each frequency. In any case, the spectral indices derived in the 
regions close to the borders of the source should not be taken 
as reliable, owing to the different sensitivities of the observing 
arrays at different frequencies and the stronger effect of noise in 
the regions of the source where the intensity is low. 

As previously reported in Bietenholz et al. (2004), the region 
around the core of M 8 1 * at each frequency has a flat (or even 
inverted) spectrum, which becomes steeper through the jet as 
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Fig. 11. Spectral-index images of M 8 1 * for different pairs of frequencies and epochs. Epoch indicates in this case the average date 
of observation (Table 1, Col. 1) of the epochs at both frequencies in each image. 



the distance to the core increases. For most of the epochs, the 
gradient of spectral index is aligned with the direction of the jet. 
This spectral-index gradient along the jet indicates how the 
opacity is decreasing with the distance to the jet origin. The same 
opacity effect indeed explains the core shift between different 
observing frequencies (reported in Sect. 14. U . In the next section, 
we discuss in detail the interpretation of our results in the frame 
of this jet model. 

5.2. M81 * and the standard relativistic jet model 

If the magnetic field in the jet and the gradient in the electron en- 
ergy distribution are large enough, synchrotron self-absorption 
can be relevant. According to the relativistic jet model, the opac- 
ity is changing along the relativistic jet and it becomes transpar- 
ent at a distance from the jet origin that depends on the observing 
frequency (Blandford & Konigl 1979). Therefore, it affects the 
apparent position of the AGN jet base (or source core) at differ- 
ent frequencies. The dependence of the core position with ob- 
serving frequency, the core-shift effect as it is known nowadays, 
was serendipitously discovered by Marcaide & Shapiro (1984) 
and later found in many sources (e.g., Kovalev et al. 20(58]). If 
the strength of the magnetic field in the jet, B, and the particle 
density, N, can both be described as power laws of the distance 
to the jet origin, the core shift, Ar, between a given pair of ob- 
serving frequencies, vi and V2, is given by (Lobanov 1998) 



Ar = Q 



vi - V2 

ViV 2 



(1) 



of the decreasing B and N. Larger Q translate into larger shifts of 
the cores at different frequencies. If there is particle-field energy 
equipartition in the jet, then k — 1. The value of Q. is the same 
for any pair of frequencies, vi and V2, as long as the power laws 
that characterize the evolution of B and N hold in the jet region 
where the cores at both frequencies are located. 

The average shift in the brightness peak between two fre- 
quencies translates into an estimate of Q.. We show in Table [5] 
the values of Q computed for each pair of frequencies shown in 
Table [3] Row 1, and assuming k = 1. All values of Q. are com- 
patible, although the uncertainties are quite large in all cases 
(of course, the uncertainties are larger for pairs of closer fre- 
quencies, since the shifts are smaller in these cases; see Eq. 
[TJ. The weighted average of Q for all the frequency pairs is 
Q= 1.75 ±0.20 mas GHz. 



5.2.1 . Position of the central AGN engine 



where Q (which could be interpreted as the core shift normalized 
to any pair of observing frequencies) depends on the power laws 



Bietenholz, Bartel, & Rupen (2000) have estimated the location 
of the central engine of M 8 1 * based on the most stable position 
of their geometrical model of M 8 1 *, referenced to the shell cen- 
ter of SN 1993J. Later, Bietenholz, Bartel, & Rupen J20041 fol- 
lowed an independent approach to determine the core position, 
based on multifrequency relative astrometry also referenced to 
the shell center of SN 1993J. Both approaches resulted in com- 
patible estimates of the location of the central engine in M 8 1 * . 
From the average value of Q. reported in the previous section, 
we can estimate the location of jet origin of M81* based on a 
larger dataset than used in Bietenholz, Bartel, & Rupen (2004), 
and using an independent approach (i.e., using Eq.Q]to find the 
jet origin). From Eq. [1] and assuming k — 1, it can be shown 
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Table 5. Q. values for the jet of M 8 1 *, derived from different pairs of frequencies. 



il (mas GHz) 



1.7-5.0 



1.7-8.4 



Frequency pair (GHz) 

5.0-8.4 1.7-2.3 



2.3-5.0 



2.3-8.4 



1.49 ±0.55 1.57 + 0.41 2.17 ±1.63 2.04 ±1.09 1.78 ±0.58 1.88 ±0.30 



Table 6. Separation between the jet origin (i.e., v — » oo) and the 
core observed at each frequency. 



Ar core (mas) 



1.7 



Frequency (GHz) 

2.3 5.0 



8.4 



1.02 ±0.10 0.76 ±0.07 0.35 ± 0.03 0.21 ± 0.02 



that the separation from the jet origin (i.e., the central engine of 
M 8 1 *) to the core observed at a given frequency v is 

Ar core = n/v, 

which translates into the separation estimates given in Table 
[6] The position of the jet origin obtained from our approach 
(Table [6]) is in excellent agreement with the position reported 
in Bietenholz, Bartel, & Rupen d2000 l l2004l . 

5.2.2. Magnetic field in the jet and black-hole mass 

Since we know the distance to M81* (3.63 Mpc; Freedman et 
al. 1 1994| l, we can compute the normalized core shift on a linear 
scale: Q ~ 0.031 pcGHz. This value can be used (together with 
the distances from the cores to the central engine of the AGN, 
reported in the previous subsection) to estimate the magnetic 
field in the jet assuming particle-field energy equipartition (see, 
e.g., Eq. 10 in Lobanov 1998). For that estimate, we also need to 
know the viewing angle of the jet (~14 degrees; Devereux et al. 
2003) and its opening angle (a projected value of ~17 degrees; 
see Sect. 154b. 

Using Eq. 10 of Lobanov (1998), the estimated magnetic 
field is 7, 10, 21, and 34 mG in the region of the brightness peak 
at 1.7, 2.3, 5.0, and 8.4 GHz, respectively. (The magnetic field is 
assumed to be inversely proportional to the distance to the AGN 
central engine.) Since these are just rough estimates, only their 
order of magnitude is considered in our discussion. 

We notice that these estimates of the magnetic field only de- 
pend on the source geometry and the differential astrometry (i.e., 
the core shift). It is intriguing that these estimates are several or- 
ders of magnitude lower than the magnetic fields derived from 
the modeling of the M 8 1 * spectrum (tens of Gauss if the in- 
verted spectrum is optically thick; see Reuter & Lesch 1996). 
Much lower magnetic fields can be matched to the spectrum if 
the emission is assumed to be produced by a mono-energetic 
population of electrons (hence to be optically thin; see Reuter 
& Lesch 1996), but such a scenario is at odds with our self- 
absorption interpretation of the observed core shift. 

Furthermore, as discussed in Lobanov (1998), the magnetic 
field at 1 pc from the central AGN engine can be used to estimate 
the mass of the black hole, by direct comparison to models of 
magnetically-driven accretion disks (Field & Rogers 1993) and 
black holes surrounded by strong magnetic fields (Kardashev 
1995). Using this approach, the mass of the black hole can be 
approximated as 

M bh ~2.7x 1O 9 BiM 



where B\ is the magnetic field at 1 pc from the AGN engine. 
With this equation, we estimate a mass of ~2 x 1O 7 M . This 
estimate is within the same order of magnitude as the value of 

" 7*i x 10 7 M reported in Devereux et al. (f2003), based on spec- 
troscopic observations of the central rotating disk of gas, and the 

~ value of 5.5j Q 6 x 10 7 M reported in Schorr Miiller et al. (1201 11 . 
based on the stellar-velocity dispersion of the bulge. As a re- 
sult, our measured core shifts, the resulting magnetic field in the 
jet, and the derived black-hole mass are all self-consistent in the 
frame of the standard jet-interaction model. 



5.3. Jet precession 

The model-fitting results shown in Fig. [2] (right) hint at position 
angles of the core Gaussians that are systematically different at 
different epochs and frequencies. From 1994 to June 2001, the 
position angle seems to increase from ~50 degrees to ~70 de- 
grees. (This increase in position angle with time seems to be 
clearer from June 1997 onwards, since the data is less spread 
during the flare from June 1997 to Jun 2001.) The change in 
position angle seems to slow down at later epochs or, even, to 
slightly reverse. 

If the major axis of the fitting Gaussian is aligned to the lo- 
cal direction of the jet in the region where the peak emission is 
produced (the most straightforward interpretation of our simpli- 
fied geometrical model), the different position angles at different 
epochs point to clear evidence of an (evolving) bent jet. In the 
following subsections, we analyze the morphological evolution 
found in M 8 1 * in more detail and discuss its interpretation as 
due to a precessing jet. 

5.3.1. Evidence of precession from the multifrequency 
astrometry 

We computed the time average of the fitted Gaussian parameters 
shown in Figs. [2] and [3] for all frequencies and for different time 
ranges. The results are given in Table [7] (the uncertainties are 
computed as the standard deviations in the averages), together 
with the average, for the same epochs, of the peak positions re- 
ported in Sect. I4.2I The quantities given in Table [7] were used 
to generate Fig. Q~2] Thus, it is a schematic representation of the 
jet structure, based on the astrometry results reported in Sect. 
I4.2I and the model-fitting results reported in Sect. [3] For clarity 
reasons, the contours shown at each frequency correspond to 0.5 
times the full width at half maximum (FWHM) of the Gaussians 
(Fig.H2(aMd)). 

A first conclusion from Fig. [12] is that, no matter the epoch 
range used in the average, the fitted Gaussians at 1.7 GHz are 
always missaligned (towards the southeast) with respect to the 
Gaussians at the other frequencies. Another conclusion is that 
the position and inclination of the average Gaussians at different 
frequencies seems to depend on the time range used in the aver- 
age. Before June 1997, the major axes of the Gaussians seem to 
be well aligned with the overall direction of the jet (Fig. 1121 a)): 
however, between June 1997 and June 2001 (the time range 
where higher peak flux densities are obtained at higher frequen- 
cies, see Sect. [3]), the Gaussians are shifted towards the east. In 
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Fig. 12. Astrometric reconstruction of the jet structure. The ellipses are contours of the Gaussians with the parameters given in Table 
|7] These parameters are the epoch-average of the model-fitting results (Table 1) for different time ranges: (a) up to June 1997; (b) 
June 1997 through June 2001; (c) from June 2001 onwards; and (d) all epochs. In (a)-(d) the contours shown correspond to 0.5 
times the FWHM. (e) is equivalent to (d), but the contours shown correspond to the FWHM. The coordinate origin in each figure 
corresponds to the centroid of the Gaussian at 5 GHz in (b). 



other words, the position angle of the overall jet increases with 
time. 

This increase in position angle with time is directly related 
to the decrease of A<S, shown in Table[7| at 1.7, 2.3, and 5.0 GHz. 
It must be noticed that the uncertainties in A<S are similar to the 
changes in Ad for the different epoch ranges used in the average, 
so the increase in the overall position angle of the M 8 1 * jet seen 
in Fig.[T2](i.e., taken from TablefT) might not be statistically sig- 
nificant. Nevertheless, we must also notice that a decrease with 
time in Ad is observed in all cases at 1.7, 2.3, and 5.0GHz (the 
increase at 8.4 GHz is only 0.16cr), thus giving more statistical 
significance to the hypothesis of an increasing position angle. 
We elaborate on this by assuming that the null hypothesis (i.e., 
a jet with a non-evolving position angle) is true. The average in 
A5 at all frequencies, Av(A£), computed from the results shown 
in Table [7] is 0.122 ± 0.056 mas. Therefore, the Student's ratio 
(see, e.g., Mandel[1964|) is 



V« Av(A<5) 



5.32, 



0~g 



where n is the number of pairs of A5 estimates (6 in our case) and 
o~s is the standard deviation of A5. The Student's ratio is the de- 



viation (in units of cr) from the average value predicted using the 
null hypothesis. Higher values of this ratio translate into lower 
probabilities for the null hypothesis. If Ad follows a Gaussian 
distribution, our set of changes in A6 follows a Student's dis- 
tribution, so the likelihood of the null hypothesis can be com- 
puted from the Student's density function using our estimate of 
the Student's ratio. The likelihood of the null hypothesis (i.e., a 
nonevolving jet orientation) turns out to be 0.08%. The likeli- 
hood of the hypothesis of a jet with an increasing position angle 
is then 99.92%. Hence, our astrometry results can be considered 
as strong evidence of an evolving position angle of the jet of 
M81*. The fitted proper motions shown in Table [4] also point, 
indeed, to an increasing position angle in the jet, since it is ori- 
ented in the northeast direction and the proper motions in dec- 
lination are negative at 1.7, 2.3, and 5.0GHz. Moreover, from 
their VSOP observations at 5 GHz, Bartel & Bietenholz (2000) 
found hints of a helical structure in the M 81* jet, which can be 
taken as a further support for a changing jet orientation. 

Therefore, we find a relation between the increasing position 
angles of the model core Gaussians (Fig. [2] right) and the in- 
creasing position angle of the overall jet, as seen from the differ- 
ential astrometry at different frequencies (i.e., Fig. [TZf aVfb) and 
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Table 7. Epoch-average of Gaussian model parameters. 



Parameter 




Epoch range 








Before June 1997 


June 1997 -June 2001 


After June 2001 


All data 






1.7 GHz 






Ao- (mas) 


- 


0.64 ±0.12 


0.68 ± 0.26 


0.65 ±0.18 


AS (mas) 


- 


-0.01 ±0.05 


-0.08 ± 0.02 


-0.04 ± 0.05 


PA (deg.) 


- 


60.8 ± 1.8 


68.8 ±3.8 


64.7 ± 4.9 


Major (mas) 


- 


2.29 ± 0.09 


2.18 + 0.11 


2.24 ±0.11 


Minor (mas) 


- 


0.94 ±0.12 


1.02 ±0.21 


0.98 ±0.18 






2.3 GHz 






Ao- (mas) 


0.45 ± 0.05 


0.43 ± 0.03 


- 


0.44 ± 0.02 


AS (mas) 


0.24 ± 0.05 


0.09 ±0.17 


- 


0.16 ±0.15 


PA (deg.) 


62.7 ± 4.0 


60.9 ± 5.3 


- 


61.5 ±4.4 


Major (mas) 


1.03 ±0.10 


1.92 ±0.07 


- 


1.62 ±0.42 


Minor (mas) 


0.5 ± 0.3 


0.5 ± 0.3 


- 


0.51 ±0.24 






5.0 GHz 






Aa (mas) 


0.09 ± 0.05 


0.09 ± 0.07 


-0.00 ± 0.28 


0.08 ±0.10 


AS (mas) 


0.02 ± 0.06 


-0.05 ±0.11 


-0.17 ±0.19 


-0.02 ±0.11 


PA (deg.) 


53 ± 12 


63.4 ± 6.2 


65.2 ± 4.5 


60.9 ±9.5 


Major (mas) 


0.82 ± 0.29 


0.93 ±0.11 


0.94 ±0.14 


0.90 ±0.19 


Minor (mas) 


0.28 ±0.15 


0.33 ± 0.06 


0.28 ± 0.03 


0.30 ±0.09 






8.4 GHz 






Ao- (mas) 


-0.06 ± 0.06 


-0.06 ± 0.07 


- 


-0.06 ± 0.06 


AS (mas) 


-0.07 ± 0.05 


-0.06 ± 0.04 


- 


-0.07 ± 0.05 


PA (deg.) 


49.5 ± 6.5 


58.5 ±5.7 


- 


52.7 ±7.6 


Major (mas) 


0.46 ± 0.08 


0.58 ±0.13 


- 


0.50 ±0.12 


Minor (mas) 


0.15 ±0.07 


0.19 ±0.04 


- 


0.17 ±0.06 



() The shifts in right ascension (Aa) and declination (A<5) refer to the coordinates given in Sect. 14.21 PA is the position angle of the core Gaussian 
(north through east direction). Major and Minor are the lengths of the Gaussian major and minor axes, respectively. 



Table |7). These results suggest a changing orientation in the jet 
of M 8 1 * (e.g., a precessing jet). There is, indeed, a possible ad- 
ditional relation between the changing orientation of the jet and 
the increase in the peak flux density observed at high frequencies 
(Fig.|2]left). A jet with an evolving orientation might help for ex- 
plaining the peak-intensity evolution shown in Fig.[2]as due, for 
instance, to changes in the relativistic Doppler boosting of the 
synchrotron emission by the plasma, as it travels through the jet, 
or to the interaction of the jet plasma with inhomogeneities in 
the surrounding medium as the jet orientation evolves. In both 
cases, a periodic evolution in the flux density of M 8 1 * would be 
expected. 

Indeed, periodic flux-density variability (with periods up to a 
few years) found in several AGN have been interpreted as related 
to jet precession (e.g., the Seyfert galaxy IIIZw2, Brunthaler et 
al. [20031 or quasar B0605 -085, Kudryavtseva et al. EOTlT) . In 
Kudryavtseva et al. (120111) . the variability (with a period 7.9 ±0.5 
years) could be related to the helical path of a jet feature. The 
period reported for B0605-085 is similar to the precession period 
that we report for M 8 1 * in Sect. I5T3~21 

It is intriguing that the flux-density flare in M 81* at 5.0 and 
8.4 GHz is not present at 1.7 GHz. However, there is a hint of 
higher flux densities at 1.7 GHz before year 2001 (values around 
0. 1 1 Jy) than after year 2001 (values around 0.08 Jy; see Fig.[2j). 
This increase in flux density at 1 .7 GHz is not as strong as the 
increases at 5.0 and 8.4 GHz (which are of about a factor 2) and 
could be indicative, for instance, of a smaller Doppler boosting 
at 1.7 GHz due to a bent jet (i.e., the viewing angle of the part 
of the jet with maximum emission at 1.7 GHz would be larger). 
Indeed, the misalignment of the core at 1 .7 GHz with respect 
to those at the higher frequencies (see Fig. [T2l might also be 
related to such a jet bending. This suggests a strong bend in the 
jet orientation within a region of ~ 0.75 - 1.0 mas from the jet 
origin (2600-3600 AU of projected linear scale). Such a bending, 



however, would conflict with the straight jet reported in Sects. 
I3.1l and [3~2l unless the bending is corrected at distances from the 
jet origin farther than for the core at 1 .7 GHz. 

Finally, in Fig.[T2Te) we show the average of the fitted mod- 
els for all epochs, but plotting the contours of the Gaussians at 
their FWHM. We can compare this figure directly to Fig. 2 of 
Bietenholz et al. ((2004). The similarity between both figures in- 
dicates similar astrometric results in both works. It must be no- 
ticed that the FWHM of the average Gaussians at 2.3, 5.0, and 
8.4 GHz are remarkably well aligned, being the location of the 
southern side of their major axes - associated with the core re- 
gion at the different frequencies - almost coincident (as reported 
in Bietenholz et al. 120041) . 



5.3.2. Evidence of precession from the evolving inclination of 
the core 



We can also study the evolution in the position angles of the 
fitted core Gaussians at each frequency. Although the spread in 
the results shown in Fig. [2] (right) is large, we notice that the re- 
sults corresponding to epochs later than year 1997 are less scat- 
tered than those at earlier epochs. This is not surprising, since 
the overall quality of the data and the sensitivity of the interfer- 
ometric arrays (i.e., amount of available radio telescopes, total 
observing time, and/or quality of the receiving systems) was im- 
proved year after year. If we focus on the subset of epochs at 1 .7 
and 5 GHz (i.e., the frequencies with the best array sensitivities) 
taken after year 1997, we end up with the results shown in Fig. 
[T3l Although we notice that there is a certain risk of overinter- 
preting these results, an oscillating trend in the position angles 
at 5 GHz seems clear. In the same figure, we show a fitting phe- 
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Table 8. Best-fit parameters for the simple precession model 
given in Eq. [2] 



Parameter 


Value 


9 (deg.) 


62.9 ± 0.3 


A (deg) 


6.7 + 0.3 


T (years) 


7.27 ± 0.08 


j5 (degyear -1 ) 


0.54 ± 0.07 


f5GHz 


1 January 1996 


fsGHz - fi.7GHz (years) 


1.9 ±0.4 



nomenological model consisting of a sinusoidal change of the 
position angle plus a small drift, i.e., 



.tMil|y(f-? v )j I- A/ /,,). 



(2) 



where V is the position angle of the core Gaussian at frequency 
v, t the observing epoch, A is the amplitude of oscillation, T 
the period, t v is the reference time at frequency v, and ji a 
small position-angle drift. We notice that the only frequency- 
dependent parameter in this simple model is t v . We give in Table 
[8] the best-fit values of the model parameters of Eq.[2] A simple 
interpretation of this model is a precessing jet of period T plus a 
frequency-dependent core shift (which accounts for t v ). We no- 
tice that Fig.[13]is compatible with the results shown in Figs. [12] 
and [10] in the frame of our jet-precession model, since the in- 
crease in position angle of the core Gaussian at 5 GHz, owing to 
the sinusoidal evolution shown in Fig. [13] begins around the end 
of year 1997 and stops around the end of year 2001. 

Thus, this increase in position angle would account for the 
rotation of the overall jet to the East (Fig. [TSl b)). If our inter- 
pretation of Fig. [13] is correct, and the long flare at 5 GHz and 
8.4 GHz shown in Fig. |2](left) is related to the jet precession, we 
would expect to see further similar flares around years 2006 and 
2015. 

Similar oscillating trends in the position angles of AGN jets 
have been found. In particular, Lobanov & Roland (2005) found 
an oscillation period of 9.5 years in the inner part of the 3C 345 
jet, which was modeled in addition to a long-term drift of 0.4 
degyear -1 . These results are very similar to the fitted parame- 
ters shown in Table [8] although the amplitude of oscillation in 
the case of 3C 345 is around 40 degrees, much larger than that 
in M81*. These authors could successfully model their set of 
22 GHz VLBI observations and multiband lightcurves with a bi- 
nary black-hole model plus a precessing accretion disk around 
the primary. The period and long-term drift in the evolution of 
the position angle of the core of M 8 1 reported here are compa- 
rable to those reported by Lobanov & Roland (2005) for 3C 345, 
but this is the first time that such behavior has been detected in 
an LLAGN. In addition, we have been able to relate the (multi- 
frequency) changing position angle of the core Gaussians to the 
(multifrequency) evolving core positions in the frame of the host 
galaxy, via the astrometry referred to SN 1993J. 



5.4. Opening angle of the jet and precession tilt 

Another conclusion that can be extracted from Fig. Q~2] is that 
the jet becomes progressively wider as the distance to the core 
increases, indicating a conic-like section in its longitudinal di- 
rection. The (projected) opening angle of the jet, estimated from 
the contours shown in Fig. [T27 a), is ~17 degrees. Since this an- 
gle has been derived from the projection of the jet on the sky 
plane, it is an upper limit to the true opening angle of the jet. 
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Fig. 13. Position angle of the core Gaussians at 1 .7 and 5 GHz, 
observed since year 1996 (typical uncertainties are given in 
Table[2]l. Predictions of our simplified precession model (see Eq. 
[2]and Table 8) are shown in continuous and dashed trace. 



If we assume that the jet is almost perpendicular to the plane 
of the rotating parsec-scale disc of gas reported in Devereux et 
al. J2003), the jet viewing angle is ~ 14 degrees. Hence, the true 
opening angle of the jet would be ~ 4 degrees. 

In the scenario of a precessing jet, the angle between the 
precession axis and the jet (i.e., the tilt angle of the jet) would be 
around a half of the total change in position angle of the major 
axes of the fitted core Gaussians (which is ~1 degrees, see the 
value of parameter A in TableO. Therefore, the real tilt angle in 
the precession (i.e., correcting for the jet viewing angle) would 
be ~1 degree. 



6. Conclusions 

We have reported on a set of global VLBI observations of the 
AGN in M81 (M81*), performed between years 1993 and 2005 
at the frequencies of 1.7, 2.3, 5.0, and 8.4 GHz. We studied the 
morphological evolution of the source and its differential astrom- 
etry with respect to SN 1993 J (a stable position reference in the 
galactic frame of M 8 1 ). 

The source consists, at all frequencies, of a slightly resolved 
core component plus a weak jet extension in the northeast di- 
rection (position angle of about 65 degrees). The core compo- 
nent has an increasing size with decreasing frequency (from 
~ 0.5 mas at 8.4 GHz to ~ 2.2 mas at 1.7 GHz). At 1.7 GHz we 
also find a more extended emission located at~ 15 mas from the 
core and at a position angle of ~ 65 degrees. This extended emis- 
sion hints at a rather straight jet in M81* at least up to 0.26 pc 
(1.1 pc) of projected (unprojected) distance to the core. The core 
components at all frequencies are elongated roughly in the di- 
rection of the jet extension and show evidence of an increasing 
position angle. 

The evolution of the peak flux density of the core shows a 
clear long-term trend at 5.0 and 8.4 GHz, with an increase start- 
ing around June 1997, a maximum reached around July 1998, 
and a decrease ending around June 2001. The flux-density vari- 
ability during this long-term flare is about 100% with respect to 
that in the quiescent stage. 
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For the astrometric study, we find a shift in the brightness 
peak of M 8 1 * towards the southwest direction with increasing 
frequency, as previously reported in Bietenholz, Bartel, & Rupen 
(2004). We compared these shifts with the predictions of the rel- 
ativistic jet model, assuming a power-law behavior for the decay 
in the magnetic field and the particle density distribution along 
the jet. We used the Lobanov ( 1998 ) model to estimate the po- 
sition of the central engine of M 81* (i.e., the origin of the jet), 
relative to the average position of the brightness peaks of our 
VLBI images. The estimated position of the jet origin is in ex- 
cellent agreement with the one reported in previous publications 
using different approaches (Bietenholz, Bartel, & Rupen 2000, 
2004). We also used the model of Lobanov ( 1998 ) to estimate the 
magnetic field in the jet and the mass of the central black hole. 
We obtained magnetic fields of a few milli-Gauss at distances 
of a fraction of a parsec from the central engine of the AGN. 
These values are much lower than those needed to fit the (self- 
absorption dominated) inverted spectrum of M81*. However, 
we obtained a black-hole mass of ~2 x 10 7 M , comparable to 
previously reported estimates using different approaches. 

From the spectral-index images of M81*, we find a region 
with inverted spectrum (typical values of a between 0.5 and 1.5) 
around the intensity peak at each frequency. The spectrum be- 
comes steeper along the jet. This behavior, together with the 
multiwavelength astrometric results, agrees with the relativistic 
jet model. 

The position of the peak intensity at 8.4 GHz is very stable 
in the frame of the host galaxy, with a proper motion consistent 
with zero to better than 10/ias year -1 . The positions of the peaks 
at the other frequencies are less stable, and they systematically 
evolve towards the south (i.e., larger position angles) at a rate of 
10, 76, and 34/^asyear -1 at 1.7, 2.3, and 5.0GHz, respectively. 
We interpret these results as strong evidence of jet precession 
towards the south. 

There is also a relation between the evolution in the rela- 
tive astrometry among the peaks at different frequencies and 
the evolving inclination of the cores, as fitted with elliptical 
Gaussians. This relation lends further support to an evolving ori- 
entation of the overall jet of M81*. The long-term variability 
found in the peak flux densities at 5.0 and 8.4 GHz also seems re- 
lated to the evolving jet orientation. In such a case, if the chang- 
ing orientation of the jet is periodic (i.e., the jet is precessing, so 
that its position angle will eventually turn northward after reach- 
ing its maximum value, around 70 degrees, and then turn back 
eastward), we predict that the flux density of M 8 1 * at the higher 
frequencies will eventually increase and follow a light curve sim- 
ilar to our observations by the next 10 years. Thus, a continued 
monitoring of the flux density and jet structure evolution of this 
LLAGN will be necessary to further confirm (and, eventually, 
better constrain) our precession model. 

An increase in the accretion rate of the central engine and 
its later correlation with the jet particle density seem unlikely to 
explain the light curve of Fig.|2](left), since it would not explain 
the simultaneous change in the jet orientation. Therefore, one 
would be tempted to think of a correlation between the flare and 
the evolving jet orientation, which would, in principle, imply the 
periodic generation of similar long-term flares in the future, as 
in the case of B0605-085 (Kudryavtseva et al. [20TTb . 
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Table 1. Results of astrometry and model fitting of M 81*. In the cases when it was not possible to perform a reliable astrometry (due to poor 
uv-coverage, too noisy data, or both), we added "-". 



Epoch 


Aa" 


A(5 a 


cr a 


Peak" 


PA b 


Major b 


Minor b 


(dd/mm/yy) 


(mas) 


(mas) 


(mas) 


(mJybeairr 1 ) 


(deg) 


(mas) 


(mas) 








1.7GMz 








06/11/93 


- 


- 


- 


0.086 


-47.0 


1.552 


0.637 


17/12/93 


- 


- 


- 


0.055 


42.0 


4.293 


0.0 


15/03/94 


- 


- 


- 


0.05 


23.7 


4.225 


0.0 


30/11/98 


0.663 


-0.052 


0.09 


0.106 


58.9 


2.324 


0.842 


06/06/99 


0.786 


-0.054 


0.09 


0.118 


59.5 


2.178 


0.919 


28/09/99 


0.598 


0.050 


0.09 


0.116 


61.3 


2.363 


1.139 


20/11/00 


0.440 


0.030 


0.10 


0.115 


63.5 


2.311 


0.858 


26/11/01 


0.741 


-0.07 


0.11 


0.073 


65.3 


2.375 


1.306 


17/11/02 


0.276 


-0.116 


0.12 


0.085 


73.6 


2.115 


1.079 


06/06/04 


0.866 


-0.105 


0.14 


0.079 


71.1 


2.085 


0.990 


06/11/05 


0.975 


-0.045 


0.15 


0.080 


64.5 


2.162 


0.712 








2.3 GHz 








19/12/95 


- 


- 


- 


0.092 


50.4 


1.512 


0.0 


13/12/96 


0.448 


0.239 


0.050 


0.103 


62.7 


1.026 


0.526 


15/11/97 


0.408 


0.214 


0.050 


0.098 


55.6 


1.913 


0.211 


07/12/98 


- 


- 


- 


- 


- 


- 


- 


13/11/00 


0.47 


-0.155 


0.080 


0.095 


66.3 


1.927 


0.797 








5.0 GHz 








06/11/93 


- 


- 


- 


0.101 


26.0 


0.787 


0.0 


17/12/93 


- 


- 


- 


0.080 


78.3 


5.224 


0.0 


28/01/94 


- 


- 


- 


0.087 


75.1 


1.541 


0.0 


15/03/94 


- 


- 


- 


0.065 


11.8 


0.021 


0.0 


22/04/94 


- 


- 


- 


0.086 


45.8 


1.148 


0.0 


22/06/94 


- 


- 


- 


0.095 


-34.0 


0.873 


0.165 


30/08/94 


- 


- 


- 


0.060 


73.2 


2.089 


0.0 


31/10/94 


- 


- 


- 


0.093 


26.7 


0.582 


0.0 


23/12/94 


- 


- 


- 


0.082 


10.7 


0.472 


0.0 


11/05/95 


- 


- 


- 


0.095 


26.0 


0.41 


0.133 


18/08/95 


0.122 


0.123 


0.023 


0.101 


44.9 


0.717 


0.0 


19/12/95 


0.126 


0.047 


0.025 


0.100 


41.1 


0.499 


0.140 


28/03/96 


0.128 


-0.016 


0.03 


0.110 


67.1 


0.759 


0.271 


08/04/96 


0.128 


-0.022 


0.03 


0.106 


61.1 


0.607 


0.220 


17/06/96 


0.088 


-0.011 


0.03 


0.104 


60.6 


0.799 


0.242 


01/09/96 


0.045 


0.001 


0.03 


0.116 


56.6 


0.901 


0.245 


22/10/96 


0.148 


0.181 


0.03 


0.096 


46.6 


0.912 


0.357 


13/12/96 


0.161 


0.073 


0.03 


0.136 


55.5 


1.215 


0.226 


25/02/97 


0.138 


0.006 


0.03 


0.105 


60.6 


1.321 


0.669 


07/06/97 


0.068 


0.021 


0.03 


0.113 


59.7 


0.855 


0.291 


15/11/97 


-0.054 


-0.103 


0.04 


0.105 


58.0 


0.829 


0.338 


18/02/98 


0.125 


0.002 


0.04 


0.123 


58.3 


0.866 


0.191 


30/05/98 


0.012 


0.006 


0.04 


0.191 


57.0 


0.827 


0.264 


03/06/98 


0.009 


0.017 


0.04 


0.154 


58.3 


0.841 


0.334 


20/11/98 


0.151 


0.032 


0.04 


0.224 


55.7 


1.084 


0.28 


23/11/98 


0.167 


0.085 


0.04 


0.229 


55.3 


1.08 


0.264 


10/06/99 


0.041 


-0.084 


0.05 


0.144 


63.5 


0.934 


0.368 


16/06/99 


0.090 


-0.400 


0.05 


0.145 


66.7 


0.891 


0.408 


22/09/99 


0.138 


-0.05 


0.05 


0.202 


69.9 


0.903 


0.386 


24/11/99 


0.249 


0.014 


0.05 


0.185 


67.8 


1.043 


0.345 


06/06/00 


0.061 


-0.140 


0.05 


0.169 


64.1 


1.128 


0.387 


24/11/00 


0.085 


-0.158 


0.05 


0.133 


69.3 


1.034 


0.326 


14/02/01 


0.167 


-0.127 


0.06 


0.106 


72.1 


0.805 


0.419 


10/06/01 


0.118 


-0.135 


0.06 


0.121 


74.8 


0.754 


0.345 


18/11/01 


0.207 


-0.124 


0.06 


0.097 


72.0 


0.867 


0.322 


24/05/22 


0.153 


-0.185 


0.06 


0.123 


72.3 


0.947 


0.229 


01/06/03 


0.463 


-0.356 


0.07 


0.121 


63.7 


0.911 


0.280 


29/03/03 


-0.203 


-0.180 


0.07 


0.068 


66.5 


1.102 


0.263 


25/05/04 


-0.103 


-0.149 


0.07 


0.098 


60.7 


0.735 


0.276 


31/10/04 


-0.046 


-0.449 


0.07 


0.105 


60.4 


0.912 


0.247 


11/06/05 


-0.498 


0.337 


0.08 


0.111 


61.1 


1.212 


0.263 


22/10/05 


-0.298 


-0.195 


0.08 


0.107 


65.4 


0.845 


0.338 








8.4 GHz 








04/08/93 


- 


- 


- 


0.101 


-71.2 


1.233 


0.655 


19/09/93 


- 


- 


- 


0.069 


-6.9 


0.291 


0.0 


06/11/93 


- 


- 


- 


0.022 


52.0 


0.337 


0.245 
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Table 1. continued. 



Epoch 


Aff a 


A<S [1 


cr a 


Peak" 


PA b 


Major" 


Minor b 


(dd/mm/yr) 


(mas) 


(mas) 


(mas) 


(mjybeairr 1 ) 


(deg) 


(mas) 


(mas) 


28/01/94 


- 


- 


- 


0.050 


37.1 


0.799 


0.0 


15/03/94 


- 


- 


- 


0.095 


39.0 


0.326 


0.103 


22/06/94 


- 


- 


- 


0.098 


51.2 


0.395 


0.142 


12/02/95 


-0.067 


-0.075 


0.02 


0.110 


56.2 


0.341 


0.185 


11/05/95 


-0.038 


-0.119 


0.03 


0.101 


49.3 


0.439 


0.056 


18/08/95 


-0.052 


0.117 


0.03 


0.098 


48.8 


0.321 


0.0 


19/12/95 


-0.194 


-0.188 


0.03 


0.116 


51.7 


0.576 


0.224 


08/04/96 


-0.039 


-0.061 


0.03 


0.110 


53.7 


0.55 


0.126 


01/09/96 


-0.137 


-0.074 


0.04 


0.112 


53.5 


0.415 


0.110 


13/12/96 


-0.176 


-0.118 


0.04 


0.133 


32.2 


0.416 


0.301 


07/06/97 


-0.120 


-0.058 


0.04 


0.109 


51.2 


0.435 


0.147 


15/11/97 


-0.113 


-0.069 


0.05 


0.158 


54.7 


0.464 


0.146 


03/06/98 


-0.035 


-0.016 


0.05 


0.205 


56.6 


0.665 


0.179 


25/02/00 


0.078 


-0.154 


0.06 


0.193 


63.3 


0.791 


0.247 


13/11/00 


-0.017 


-0.016 


0.06 


0.146 


66.7 


0.54 


0.223 



a Aa and AS are the right ascension and declination (respectively) of the intensity peak referred to the coordinates given in Sect. 14.21 a is the 
uncertainty in both, Aa and A<5. 

b Peak, PA, Major, and Minor are the intensity peak, position angle, major axis, and minor axis (respectively) of Gaussians fitted to the M 81* 
visibilities. 
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